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'  TIME-FREQUENCY  ANALYSIS  OF  ACOUSTIC  IMAGERY 

FOR  BOTTOM  MAPPING 

Maria  T.  Kalcic*  and  Andrew  B.  Martinez^  and  Douglas  N.  Lambert^ 


Abstract — Acoustic  imagery  generated  by  &e 
Naval  Researdi  Laboratory’s  ASCS  (Acoustic  Sed> 
iment  Classification  System)  is  analyzed  using  time- 
frequency  methods.  Hie  ASCS  imagery  is  gener¬ 
ated  using  normal-incidence  acoustic  data.  Nonlin- 
earities  in  the  propagating  medium  produce  multi¬ 
ple  harmonics  in  tiie  return.  Different  components 
of  the  bottom  and  sub-bottom  produce  different 
spectral  signatures,  based  on  the  harmonics  that 
be  used  to  classify  tiie  bottom.  Although  the 
exact  acoustic  properties  of  the  media  may  not  be 
known,  it  is  possible  to  segment  the  spe<frally  dif- 
ferentcomponents.  Hiis  paper  presents  an  ef^ent 
method  for  performing  the  time-frequency  analy- 
tis.  Ihe  tradeoffs  between  time  and  freqnracy  res¬ 
olution  are  examined  and  exploited  to  produce  the 
best  overall  resolution.  Examples  using  tiie  ASCS 
data  demon^rate  Its  use  for  classification  of  acous¬ 
tic  imagery.  Buried  pipes,  methane  gas  deposits, 
and  other  features  can  be  differentiated  using  this 
method. 

I.  Introduction 

A  need  exists  for  detailed  mapping  of  the 
seafloor  prior  to  dredging  operations  and  the 
laying  of  pipelines  and  cables.  In  addition, 
there  is  a  need  to  be  able  to  precisely  locate 
existing  buried  pipeline  and  cable  assets  in 
the  area  so  that  they  are  not  damaged  during 
flie  new  operations.  The  location  and  positive 
identification  of  these  man-niade  buried  ob¬ 
jects  can  often  be  difficult  because  they  can 
be  masked  in  the  subbottom  imagery  by  natu¬ 
ral  sediment  structures,  methane  gas  pockets, 
and  man-made  discarded  trash.  However,  we 
have  found  that  the  use  of  a  high  frequency, 
short  pulse  and  nonlinearities  in  the  trans¬ 
mitting  medium  produce  a  broadband  signal 
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fliat  can  be  exploited  to  better  distinguish  be¬ 
tween  these  objects.  This  paper  describes 
how  these  properties  can  be  enhanced  using 
time-frequency  approaches. 

n.  Background 

Ibe  narrow  beam,  short  pulse  and  wide 
bandwidfli  of  the  AS(ZS  have  shown  excellent 
results  in  the  determination  of  detailed  sub¬ 
bottom  structure  and  prediction  of  geoacous- 
tic  properties  [1-3]  as  well  as  the  detection  of 
buried  objects  such  as  pipelines  in  the  upper 
few  meters  of  the  seafloor.  Fig.  1  shows  an 
oil  pipeline,  0.51  m  in  diameter,  buried  in  one 
meter  of  sand,  detected  with  the  ASCS.  The 
pipeline  has  an  extended  reverberation  tail 
due  to  ringing  inside  the  pipe.  Also  visible 
in  the  image  are  sediment  disturbances  due 
to  trenching,  seen  directly  above  the  pipeline 
in  Fig.  1  as  sort  of  ’V’  shaped  dark  area. 
Remnants  of  the  trenching  berms  are  evident 
above  and  to  either  side  of  the  trench.  Many 
of  the  returns  from  muddy  sediments,  how¬ 
ever,  are  embedded  with  methane  gas  bub¬ 
bles  that  display  a  strong  acoustic  impedance 
similar  to  those  of  the  buried  objects  caus¬ 
ing  false  alarms  if  the  system  were  used  for 
buried  object  detection.  It  is  necessary  to  de¬ 
velop  techniques  that  can  distinguish  hard  re¬ 
flectors,  such  as  gas  pockets,  from  pipelines 
and  cables,  in  order  to  automate  the  process 
of  object  location. 

ni.  Broadband  ASCS  Acoustic  Imagery 

Acoustic  imagery  of  a  sediment  cross- 
section  is  produced  with  a  single,  normal 
incident,  narrow  beam  (12  degrees)  pinging 
along  a  fixed  line.  The  ASCS  transducer  used 
in  this  study,  an  ELAC  LSE-179,  produces  a 
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F>g.  1 ,  Oil  pipeline  deiecied  by  ASCS  ai  1 5kHx  show¬ 
ing  extended  reverberation  tail,  and  sediments  dis¬ 
turbed  by  trenching  operations. 

short  1 5  kHz  pulse  that  results  in  a  broadband 
signal,  with  high  lime  resolution.  Sixteen  bit 
data  provide  approxi  mately  90  dB  of  dynamic 
range. 

The  acoustic  imagery  generated  by  the 
ASCS  contains  reflections  from  many  objects 
with  the  same  amplitude  making  them  in¬ 
distinguishable.  Often  the  problem  of  dis¬ 
tinguishing  methane  gas  deposits  from  oil 
pipelines  or  other  objects  is  very  difficult  to 
the  trained  viewer,  and  much  more  sr)  to  an 
automated  detection  system.  Fig.  2  shows  a 
log-intensity  image  of  several  hundred  pings 
stacked  as  columns  to  produce  a  single  im¬ 
age  witli  depth  increasing  from  top  to  bot¬ 
tom.  Methane  gas  pockets,  surface  debris 
and  a  buried  cable  are  visible  in  the  image. 
The  log-intensity  is  used  so  that  some  of  the 
less  reflective  materials  are  more  visible,  'fhe 
bright  object  in  die  middle  left  of  the  image 
is  a  buried  cable.  As  one  can  see,  other  ob¬ 
jects  in  the  image  appear  very  similar  and 
it  is  difficult  to  differentiate  the  cable  from 
the  methane  at  the  bottom  of  the  image  or 
seashell  debris  at  the  surfitce. 

The  shortcoming  of  this  approach  is  evi¬ 


Fig.  2.  Several  hundred  pings  from  the  ASCS  witii 
intensity  in  dB 


dent  by  looking  at  one  column  (ping)  from 
Fig.  2.  Figure  3  shows  the  broadband  in¬ 
stantaneous  intensify  of  the  time  sequence  of 
one  ping.  The  largest  peak  corresponds  to 
the  buried  cable.  The  initial  large  peak  rep¬ 
resents  seashell  clutter  at  the  surface  and  the 
peaks  near  the  bottom  (around  samples  550- 
600)  represent  methane  gas  pockets.  Since 
this  shows  the  time  varying  intensity  for  all 
frequencies,  differentiation  of  objects  based 
on  different  frequency  characteristics  is  not 
possible,  and  all  processing  must  be  purely 
spatial. 

IV.  Narrow-band  Fh'ocessing 

The  transmis.sion  of  a  31.25  kHz  signal 
with  a  15-kHz  piezometer  transducer  pro¬ 
duces  hamionics  in  the  received  signal  due 
to  nonlineaiities  in  the  transmission  medium. 
This,  combined  with  the  use  of  a  short  pulse 
(0.13  ms),  re,suUs  in  a  broadband  .spectrum 
with  .strong  peaks  at  several  frequencies.  It 
turns  out  that  difTerent  materials  respond 
quite  differently  to  the  transmitted  signal, 
suggesting  an  investigation  into  the  spectral 
character  of  the  returns. 

Fig.  4  show.s  the  spectrum  of  one  ping 


Fig.  3.  Instantaneous  broadband  intensity  of  a  single 

ping  through  area  of  interest. 

(column)  of  Fig.  2  containing  the  cable  and 
metihane  gas  pockets  in  cross-section.  The 
strong  peaks  at  various  frequencies  indicate 
that  more  information  might  be  available  by 
examining  those  frequencies  separately. 

The  approach  used  to  examine  selected 
center  frajuencies  is  to  slice  the  frequency 
spectrum  into  narrow  bandwidths  by  apply¬ 
ing  a  narrow-band  filter,  in  this  case  a  nth- 
order  forward-backward  Butterworth  filter. 
The  filter  is  designed  for  zero  phase  shift  and 
centered  at  the  frequency  of  interest.  Fil¬ 
tering  the  original  pings  and  plotting  the  in¬ 
stantaneous  intensity,  results  in  profiles  sig¬ 
nificantly  different  from  that  of  the  original 
data.  The  results  for  frequencies  at  12.5, 19.4, 
33.75,  35.0,  and  55.6  kHz  are  shown  in  Fig. 
5. 

The  profiles  show  a  strong  return  of 
methane  gas  at  low  frequencies  with  an  in¬ 
creasing  return  from  the  buried  cable  as  the 
frequencies  increase.  In  the  12.5  kHz  pro¬ 
file,  the  methane  gas  peaks  above  the  other 
returns,  with  an  intensity  of  over  100.  In  the 
1 9.4  kHz  profile  the  methane  gas  continues  to 
peak  over  100  but  the  buried  cable  begins  to 
;  pe^‘  cibse  to  100  also.  There  is  also  a  shift  in 


Fig.  4.  Power  i^jectrum  of  one  ping  through  cable  and 
methane  gas  pockets. 


the  location  of  the  methane  gas  responding  at 
this  frequency.  At  33.75  kHz  the  cable  peaks 
above  5CK)  and  the  relative  return  from  the 
methane  gas  is  strongly  reduced.  The  same 
occurs  at  35  kHz.  At  55.6  kHz  the  methane 
gas  is  no  longer  visible  relative  to  the  cable 
due  to  signal  attenuation  and  the  surface  de¬ 
bris  is  the  most  visible  due  to  surface  refla¬ 
tion  at  the  higher  frequency. 

The  application  of  the  filter  to  the  entire 
image  produces  a  set  of  images  at  different 
frequencies,  shown  for  the  area  surrounding 
the  cable  cross-section  in  Fig.  6.  These  im¬ 
ages  can  be  considered  a  multispectral  data 
set.  The  cable  cross-section  is  easily  seen 
in  the  33-35  kHz  frequency  range  and  the 
methane  is  most  visible  at  the  19.4  kHz  slice. 

There  are  several  important  facts  to  point 
out  in  these  figures,  one  being  that  the  peak 
heights  are  relative,  not  on  the  same  scale, 
due  to  the  attenuation  of  higher  frequencies 
by  the  sediment  medium.  Also,  because  these 
heights  are  relative  along  the  time  axis,  the 
careful  selection  of  time  and  scale  can  be 
combined  to  enhance  features  in  time  as  well 
as  frequency.  For  example,  by  looking  at 
a  time  slice  of  the  original  data  around  the 


;Fig.  5.  Iit$£aot9ftep:Os  jjarrpw-hand^^j  yersus 

tiiejcrs  i>elow  tbo  seabed  of  a  single  ping  a(  12.5* 
and  55:d  kHz 


methane  peak  at  55.6  kHz,  the  methane  cafi 
be  detected  as  it:h^  a  stronger  refuni  than  the 
local  sediment,  between  3  and  3.5  m.  Thirdly, 
the  diametef  of  the  cable  makes  it  almost  in¬ 
visible  at  the  lower  freqnencie.s,  .so  the  knpwlr 
edge  of  the  object’s  size  can  be  used;  to  sciect 
the  scale  or  ffeqoeticy  that  best  locates  it.  The 
combinatibn  qf  time  and  frequency  together 
boo.sts  the  detection  and  identification  power 
of  an  objection  recognition  system, 

V.  Time-Frequency  Analysis;  with  Wavelets 

Wavelets  are  a  function  of  scale  and  time, 
where  scale  is  closely  related  to  frequency. 
Two-dimen.s{onal  wavelet:  analysis  is  used 


Fig.  6,  Narrow-band  image  slice.s  at  selected  frequent 

cies; 

here  to  decompose  the  original  image  into 
su  bi  mages  where  each  pi  xel  value  represents 
the  conuibutiOn  of  it’s  lOGaliiy  to  a  particular 
scale  (or  frequency)  in  (ime.  tn  the  im¬ 
age  is  decomposed  into  horizontal,  diagonal 
and  vertical  components  at  each  level.  Each 
level  gives  a  low-p^S  component  (approxi¬ 
mation)  and  a  high-pass  component  (detail). 
The  level- 1  approximatioivis  decimated  by  2 
and  filtered  again  to  produce  a  level-2  set  of 
high-  and  jow-pass  images. 

The  ASCS  acou.stic  image  should  be  ideal 
for  wavelet  analysis.  As  demonstrated  above, 
different  components  in  the  original  image 
have  different  vertical  .spectral  lesponse  re¬ 
lated  to  their  Compb.sitioh  and  scale;  The  hor¬ 
izontal  component  of  the  image  is  strictly  a 
function  of  object  scale.  W'avelets,  by  re¬ 
taining  both  scale  and  frequency  information 
.should  tran.sform,  the  image  in  a  manner  con¬ 
sistent  with  our  concept  of  both  frequency 
and  scale;  By  selectively  thre,shordirig  co¬ 
efficients  unrelated  to  the  object  in  scale  or 
frequency,  we  Can  effectively  “de-noise”  the 
data. 

As  an  illustration  of  this,  a  fourth-order 
paubechies  wavelet  transform  was  applied  to 


Rg.  7.  ITie  level  2  horizontal  deKiil  coeflicients  of  the 
fourth  order  Daubechies  wavelet  transform  and  the 
reconstructed  image  after  de-noising. 

tjie  original  utifiltei^  im  2.  The 

results  were  quite  good.  Fig,  7  shows  {he 
levelrl  horizontai  detail  of  the  original  M* 
age,  selected  because  it  shows  the  buried  ca¬ 
ble:  very  deafly.  The  wavelet  coelTicients  may 
also  be  used  to  de-irolse  the  image  through: 
thfesholdiijg  of  tl)e  coefficients  for  each  level 
and  direction,  and  reconstfuction  of  the  im¬ 
age  from  the  Goetlicients  below  the  threshold. 
An  example  of  a:  de-noised  image  is  shown  in: 
Fig.  7, 

VI,  Conclusions 

In  this  paper,  we  have  shown  that  buried 
objects  can.bo  detected  using  both  multispeC- 
tral  and  time-:frequency  approaches. 

A  frequency  based  approach  works  well  in 
this  environment  where  the  objects  respond 
differently  to  different  frequencies.  By  slic¬ 
ing:  the  Spectrum  into  narrow  bands,  it  is 
possible::  to  distinguish  objects  by  exploiting 
these  differences.  This  frequency  slicing  a|H 
proach  h:as  the  further  advantage  of  produc¬ 
ing  a  multispectfal  data  set  which  lends  itself 
well  to  multivariate  classification  techniques. 
..Results  thus  far  have  shown  that  i.solat- 


ing:  the  object  in  frequency  and  then  in  time 
and  scale  using  wavelet  analysis  is  ah  effec¬ 
tive  detection  tool.  Knowledge  of  the  ob- 
jeef  s  acoustic  impedanesa  properties  m  well 
it’s  sl?e  can  be  used  to  define  the  .chafactef* 
istics  of  the  wavelefs  order  and  level  of  de¬ 
composition  to  automatically  look  for  objects 
m:  certain  scales. 
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